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Non-Abelian quantum holonomy of hydrogen-like atoms 
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We study the Uhlmann holonomy [Rep. Math. Phys. 24, 229 (1986)] of quantum states for 
hydrogen-hke atoms where the intrinsic spin and orbital angular momentum are coupled by the 
spin-orbit interaction and subject to a slowly varying magnetic field. We show that the holonomy 
for the orbital angular momentum and spin subsystems is non-Abelian, while the holonomy of the 
whole system is Abelian. Quantum entanglement in the states of the whole system is crucially 
related to the non-Abelian gauge structure of the subsystems. We analyze the phase of the Wilson 
loop variable associated with the Uhlmann holonomy, and find a relation between the phase of the 
whole system with corresponding marginal phases. Based on the result for the model system we 
provide evidence that the phase of the Wilson loop variable and the mixed-state geometric phase 
[E. Sjoqvist et al. Phys. Rev. Lett. 85, 2845 (2000)] are in general inequivalent. 

PACS numbers: 03.65.Vf, 03.65.Ud, 31.15.aj 



I. INTRODUCTION 

The pioneering work by Berry [l| and Wilczek and Zee 
, have triggered considerable interest in effective gauge 
structures in the adiabatic evolution of non-relativistic 
quantum systems. Non-Abelian quantum holonomies 
have been examined in the context of nuclear rotations of 
diatoms nuclear quadrupole resonance 0, semicon- 
ductor heterostructures Q, trapped atoms Q, quantum 
optics and superconducting systems [1]. It has been 
pointed out @ that non-Abelian holonomy may be used 
in the construction of universal sets of quantum gates for 
the purpose to achieve fault tolerant quantum computa- 
tion. 

In Refs. non-Abelian holonomies are related to 

the existence of degenerate energy eigenstates that can 
be controlled by a set of slowly changing parameters. 
In contrast, Uhlmann [lo| has shown that non-Abelian 
gauge structures may appear along sequences of density 
operators representing mixtures of quantum states, ir- 
respective of the degeneracy structure of the underlying 
Hamiltonian. Such non-Abelian structures may arise for 
subsystems of composite systems undergoing adiabatic 
evolution, since the marginal states are mixed if the in- 
stantaneous energy eigenstates of the whole system are 
entangled. 

The purpose of the present paper is to follow the 
Uhlmann approach to examine non-Abelian gauge struc- 
ture in the case of spin-orbit {LS) coupled hydrogcn-like 
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atoms subject to a slowly varying magnetic field. We 
show that the adiabatic Uhlmann holonomies for the spin 
{S) and orbital (L) parts become non-Abelian although 
the one of the whole LS state is Abelian. Studies of 
the same model were carried out previously [ll| (see also 
[l^ ) by using the mixed-state geometric phase approach 
developed in Ref. [ij]. In particular, in Ref. [Hi] it was 
shown that the mixed-state geometric phases of the L 
and 5* subsystems always sum up to the standard pure 
state geometric phase of the whole system. In contrast, 
we show here that the phases of Wilson loop variables as- 
sociated with the Uhlmann holonomies satisfy this sum 
rule only for specific paths, while for other paths there 
is a deviation of tt from the sum rule. This deviation 
from the sum rule demonstrates a strikin g no n-trivial 
difference between the Uhlmann holonomy [l^j and the 
mixed-state geometric phase [l3|. 

The outline of the paper is as follows. In Sec. |ll] the 
Uhlmann holonomy is briefly reviewed. The correspond- 
ing parallel transport equations for adiabatic rotation of 
angular momentum states are derived in Sec. IIIII In 
Sec. IIVI we compute the Uhlmann holonomies for the 
L and S subsystems as well as that of the total angular 
momentum. We examine in particular the non-Abelian 
nature of the subsystem holonomies as well as the ad- 
ditivity of the phases of the Wilson loop variables asso- 
ciated with the Uhlmann holonomies. The paper ends 
with the conclusions. 



II. UHLMANN HOLONOMY 

In this Section we summarize the main definitions and 
properties of the Uhlmann holonomy [l^]. Let C : [0, 1] 3 
t pt he a, smooth path of density operators acting on 
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some Hilbert space H. An operator Wt such that pt = 
WtWf is cahed an amphtude of pt- Wt can be written 



as Wt — y/ptVt, where the "phase factor" Vt is a partial 
isometry on H. For any choice of Wq = Wq, there is 
a difFerentiable path [0,1] 3 t Wt of ampUtudes over 
C that satisfy the parallel transport condition 



WUW = dW^W, 



(1) 



where dW = dtWt- Inserting W — into Eq. ([T]) 

yields [3 



dvv'^p + prfyyt = d^^ - Vp^Vp. 



(2) 



By solving for Y ^ we define the Uhlmann holonomy as- 
sociated with the path C to be 



(3) 



The operator C/uhi(C) is a unique partial isometry 
(unique unitary if all pt are full rank), gauge in- 
variant (i.e., independent of the choice of Wq), and 
reparametrization invariant (i.e., independent of the 
speed of the evolution along C); thus, C/uhi(C') is a prop- 
erty of the path C. 

In Uhlmann's approach, a system in a mixed state is 
thought to be a subsystem embedded in a larger quan- 
tum system which is in a pure state. The pure state 
is referred to as a purification of the mixed state. This 
is accomplished by introducing an auxiliary system with 
which the original system is entangled. The purification 
is equivalent to the amplitude Wt of pt- If Wt satisfies 
the condition in Eq. ^ along the path C, then, inspired 
by the pure state geometric phase [l6| . one may assign 
the phase 



(^uhi = arg Tr(WoVi) 



(4) 



to the Uhlmann holonomy J7uhi(C'). The definition (/^uhi 
for Uhlmann phase has been used to investigate theoret- 
ically [l7l-[l9| and experimentally [l^ apossible relation- 
ship between the Uhlmann holonomy [10"| and the mixed- 
state geometric phase /? [13|]. The latter differs consider- 
ably from (^uhi in that it is, for cyclic unitary evolution, 
the sum of geometric phase factors of the eigenstates 
of the density operator weighted by the corresponding 
eigenvalues. Thus, contrary to lySuhi, /? is essentially a 
decomposition dependent and Abelian geometric phase 
concept. 

On the other hand, as for the case of the non- Abelian 
Wilczek-Zee phase factor Q, the Uhlmann holonomy in 
Eq. ([3]) for cyclic evolutions takes the form of a Wilson 
loop Pe~*^c"* for a vector potential A, where P stands 
for path ordering. On the basis of this fact and on the fact 
that the Wilczek-Zee phase factor is a natural extension 
of the Berry phase [1] to systems with degenerate spectra, 
one can argue that the phase of the Wilson loop variable 



Tr ^Pe ^ic^ is more natural quantity than i^Suhi, and 
define the Uhlmann phase as 



7 = arg [Tr(C/uhi(C))] 



(5) 



Note that the two phase quantities (/^uhi and 7 need not 
be equal. The focus of this paper is to turn our attention 
to the Uhlmann phase defined in Eq. ((SJ. In particular, 
we demonstrate that 7 behaves very differently from the 
mixed-state geometric phase [l^ 



III. UHLMANN HOLONOMY OF ANGULAR 
MOMENTA 

As an introduction to the model studied in Sec. IV, 
we consider the case of adiabatic transport of a quan- 
tum angular momentum S. This angular momentum is 
assumed to be coupled to another quantum angular mo- 
mentum S^''-'. The coupling is assumed to be spherically 
symmetric. Both S and S^*"^ are exposed to an external 
classical magnetic field B. The Hamiltonian of the sys- 



tem takes the form H = IJ^' 



tott / 



where 



U'°\e,4>) = e~'^^'^' e^^'^^T e'^sf' and S*°' = S + S^"") = 
(5*°S 5'^°*, 5*°') {h = 1 from now on). 9, are spherical 
polar angles parametrizing the direction n = B/|B| of 
the external magnetic field. 

Hz is independent of the polar angles [9, (f)) of the two- 
dimensional parameter sphere of all possible magnetic 
field directions. We assume that [Sl°*,Hz] = to make 
sure that H is well defined on S^. The energy eigenstates 
can be represented by the smooth vector- valued functions 
|^(")(6l,0)) = C/t°*(6l,0)|4"^), well defined on the open 
patch 52-{6' = tt}, and IV'^"^^,^))' = f/*°*(6', 0)|4"^) = 
[/*°'(0,0)e"2«0Sf'|^(")^^ ^^l]^ defined on the open patch 

52 - {6* = 0}. Here, = i^^^^i"^). The vec- 

tors </>)) and IV''"'' (^' '/'))' define two monopole 

sections |2l| over the parameter sphere. These sections 
are related by a single-valued gauge transformation so 
that 

|V(")(0,0))(V^(") (0,0)1 = |V("H0,0))' '(^(")(^,</')l (6) 

in any overlapping region on the parameter sphere. 

The reduced density operator p(")(0, 0) representing 
the marginal state of S, corresponding to the nth energy 
eigenstate of H, is obtained by partial trace Tr^ over the 
degrees of freedom associated with S^'^\ i.e.. 



p(")(0,0) = U{9,cf,)pi"^U\9,cb) 
= Ui9,^)pi"WH9,<ly). 



(7) 



Here, U{9, 0) 



-i(pSz 



e ^^Sy^t4>Sz rotation opera- 



tor, U{9,cb) = [/(0,0)e-2^^^^ and pi"^ = Tr,|vi"^) (4"' I 
is a "reference" state, p'"' {9, (f>) defines the mixed state 
of our subsystem, for the nth energy eigenstate. 
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In the adiabatic regime, the path F : [0, 1] 9 t — > 
{6t, (j)t) on the parameter sphere 5^ of magnetic field di- 
rections maps to the path C^") : [0,1] 3 t ^ P^'^KOt^ 4>t) 
in state space of the considered angular momentum. Let 
= U{9,<j))V'''^^ be a partial isometry that satisfies 
parallel transport along F. With d — d9dg + ^03^, we 
have 

= -2iy^[d0(l - cos 61)5, 
+d(f> sm9{Sx cos (j) + Sy sin ( 

^ (8) 



-d9{—Sx sin (f> + Sy cos 4>)] \/ pi"^ 




along F. 

Repeating the calculation for the other monopole sec- 
tion, by using the decomposition V^'"' ~ U{9,(I))V'^^\ 
leads to the equation 



jT7(")t7(")^ (") I (") jT7(")T7(")t 

dV V Pz + Pz dV V 



= 2i\J pf' [d(i)(l+ cos 9) Sz 
-dip sm9{Sx cos (j) — Sy sin 



+d9{Sx sin (p + Sy cos 




along F. Since 5,] = 0, 



Tin) 



(9) 



choice Vq'^ = g2»,^oS,y^l») satisfies Eq. El Thus, the 

difference between v''"'^ and V^^^ precisely compensates 
the difference between the rotation operators U {9, (f) and 
U{9, (j)) so that the Uhlmann holonomy remains the same 
in the two representations. In other words, J7uhi(C^"^) 
is independent of which monopole section we use. This 
implies that either of the above pair of monopole sections 
can be used to calculate the Uhlmann holonomy for any 
path on the parameter sphere. 



IV. UHLMANN HOLONOMY OF 
HYDROGEN-LIKE ATOMS 



In Ref. the adiabatic geometric phases of the 

L^-coupled hydrogen atom in a slowly rotating mag- 
netic field B Bn — B {sin 9 cos cp, sin 9 sin (f), cos 9) were 
analyzed. The adiabatic geometric phases of the whole 
system and of the orbital (L) and spin (S) angular mo- 
mentum subsystems were computed. In particular, it 
was demonstrated that the subsystem phases add up to 
the phase of the whole system. The purpose here is to 
compute the corresponding Uhlmann holonomies and to 
examine their relation. 



A. Model system 



We consider the spin-orbit (LS) part 



iJ„ = gn-(L + 2S) + 2L • S 

= Uji9,(f>)HzUl{9,^), 



(10) 



of hydrogen-like atoms exposed to an external magnetic 
field pointing in a direction defined by the unit vector 
n. The first and second terms are the Zeeman and LS- 
coupling contributions, respectively, g is the Zeeman-L5 
strength ratio (assumed to be time-independent), and n 
defines the adiabatic parameter sphere under slow 
changes in the direction of the external magnetic field. 
We may choose 



Ux{9,^) 



X^L,S,J, (11) 



where J = L + S is the total angular momentum and Hz 
is the Hamiltonian at the north pole n = (0,0, 1) of the 
parameter sphere. 

The Hamiltonian Hz is block-diagonalizable in one- 
and two-dimensional blocks with respect to the product 
basis with elements |Z,m)|i,±^) = \l,m)\±) being the 
common eigenvectors of L^, L,, S^, 5,. Each block may 
be labeled by the eigenvalue p^ — + + ^ 

of Jz- The two extremal subspaces characterized by 
\p\ = 1 + -^ = Pc are one-dimensional corresponding to the 

two product vectors |i/'± '^'"')) = |;,±0|±). The remain- 
ing blocks are two-dimensional, each of which spanned by 
the vectors \l, m = p—-^)\+), \l, m = p+^)\—), \p\ < 1+-^- 

For each fx, the corresponding energy eigenvectors \4'±'^^) 
become LS entangled. The amount of LS entanglement 
may be measured in terms of concurrence 



c(v.^'^') 



where 



cos a 



2fi 



(12) 



(13) 



Agn + (21 + 1) 



Note that C*^''^^ is independent of the ± label and varies 
between for product states (a^''^) — 0, tt) and 1 for 
maximally entangled states (a^'^^-* = 7r/2). 

The instantaneous energy eigenvectors of H^ are re- 
lated to the eigenvectors of Hz according to 



= UL{9,q^)Us{9,cPM^l-^^} 



(14) 



Let F : [0, 1] 9 t — {9t, pt) be a parametrized path on 
the parameter sphere. Assume that the external mag- 
netic field slowly traverses F so that the adiabatic ap- 
proximation is valid. With this assumption F maps to 
the paths 



: [0, 1] 9 t ^ p^;,'' 
= Ux{9tAt)p^kiu],{9u 



t) 



t), X = L,S,J, (15) 
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in the spaces of density operators, where 



Pj,± 
Pl,± 
Ps,± 



Tr^l^, 



Cm) 



(16) 



In the foUowing subsections we compute the Uhlmann 
holonomies of the total angular momentum J of the atom 
and its subsystems L and S under the assumption .g 7^ 0. 



B. Holonomy of total angular momentum 



Let Vj'^-* denote the solution of Eq. ([5]) with reference 
state |V'±'^'')(V'±''^'' I- We obtain the partial isometry 

= e-'^i'r(i-™^«)'^"*|V'£'^))(V'£^^Vi.S> (17) 
which yields the Uhlmann holonomy 



g — i/^ /p(l — cos6 



(18) 



To compare with the corresponding geometric phase 

factor e*^^'^-'.± \ we note that while this geomet- 
ric phase factor is a unit modulus complex number, 
Uuh\{C^j'±) in Eq. (IT^ is a partial isometry with 

uUc'i^yUu^ic'j^) and UMCfX^)uUcfX^) be- 
ing projection operators onto the initial and final states, 
respectively. In particular, this shows that while the geo- 
metric phase /3(Cj^^) is undefined if the two end points 



correspond to orthogonal states, ?7uhi(Cj^ ) is 



of C, 

a well-defined partial isometry. On the other hand, a 
direct calculation yields 



Tr[C/,hi(C,yf )] = \TT[UMCf^^)]\e 



(19) 



that demonstrates an explicit relation between the Wil- 
son loop variable associated with the Uhlmann holonomy 
and the geometric phase factor, unless Tr[[/uhi(Cj 
vanishes, which happens precisely when the initial and 
final states are orthogonal. These results establish a one- 
to-one relation between the standard pure state geomet- 
ric phase (l6| and the corresponding Uhlmann holonomy 
of the whole system. 



C. Holonomy of the L and 5* subsystems 

Now we compute the Uhlmann holonomies of the L 
and S subsystems. Let us start with the extremal states 
fi = ±fJ.c- We note that the eigenvectors 



,0) = UL{0,m,±l)Us{O,(f')\±) (20) 



of Hn are tensor products of states of the two subsystems 
L and 5*. We thus find the Uhlmann holonomies for the 
L and S subsystems as 



[/uhi(4'±'"^) 



COS 0)d(f) 

xUL{ei,(t>i)\i,±i){i,±i\uli0o,M, 

xC/s(0i,0i)|±)(±|C/l(0o,</'o). (21) 

Note that the associated holonomies are ^-independent 
and satisfy the product relation 

c/uhKc^/i ^°^) = [/uhi(ci';^^'=)) ® uMc^si""^)- (22) 

Next, we compute the Uhlmann holonomy in adiabatic 
evolution of non-extremal energy eigenstates character- 
ized by < I + ■^^ The marginal states are rank- two 
density operators obtained by adiabatically evolving the 
states p^l''^ and Pg'^ under unitaries Ul and Us, re- 
spectively. One may solve Eq. ([5]) with reference states 

P^l'± for a path F on the parameter sphere to obtain the 
Uhlmann holonomy of the L subsystem. 



x;7i(0o,0o)- 



(23) 



Here, we have introduced the vector potential compo- 
nents 

1 „^ / 



A 



2 \ le 







d9. 



A^i:^^J = ^(l-cos0) f J \ 



(24) 



1 / — 1 + cos sin 9e^' , , , 



Wt 



smf 



The vector potential Af^^ = A 



1 — COS 

- AlI^^^ is expressed 



in the basis {|/^— ^) , and w ~ \J + \Y ^ P-"^ 

0. 



> 



It is worth noting that the vector potential ex- 
hibits a U(l) part being proportional to the identity. In 
the case where F is a loop, this part gives rise to the 



global geometric phase factor e~*'^ 
r2 being the solid angle enclosed by F on the parameter 
sphere 5^. 

Similarly, for the S subsystem we have 



C/uhi(4':± ) = t/s(^i,<^i)Pe"'^^(^ 
xC/^(0o,'/'o) 
with vector potential components 



(25) 



A 



S;e 



1 



A 



le 
2' V -ie*"* 

If I-COS6' Cf'^'^^sinf 



-i4> 



C^'^) sinefe*"^ 



-1 + COS ( 



d(j) 



(26) 
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in the {|+), |— )} basis. Note that A^g'^^ does not have a 
U{1) part. 

Unhke the extremal case, the marginal Uhlmann 
holonomies are g-dependent via the concmrence C^'''^\ 
Furthermore, there is a dimensional mismatch between 
the rank-one holonomy of the J system and the rank- 
two holonomies of the L and S subsystems; in general 
this mismatch implies that there is no path for which a 
product rule similar to that in Eq. (|22|) holds. 

It is noticed that when ^ 0,7r, i.e., when the 

concurrence C^''^^ is non-zero, the vector potentials con- 
tain non-Abelian components. In other words, the non- 
Abelian nature of the subsystem holonomies is due to 
entanglement. This is analogous to the Levay geometric 
phase defined for two-qubit systems p3j , which is a path- 
dependent unit quaternion that may find realization in 
two-particle interferometry . The holonomy group as- 
sociated with this geometric phase becomes, just as the 
Uhlmann holonomy of the subsystems, Abelian in the 
product state case. 

In the "classical" limit characterized by Z — >■ cx) and 
1^1 = 0(0, C^'-'') vanishes as a^''^) tends to zero. Thus, 
the L and 5* holonomies turn Abelian in this limit. We 
further see that for cyclic evolution in the l/\^\,l/\g\ — 
oo limit, C^''''^ tends to its maximum value 1, i.e., the en- 
ergy eigenvector becomes maximally entangled, and the 
holonomies turn into the transpose of the Wilczek-Zee 
holonomy for nuclear quadrupole resonance setup dis- 
cussed in Ref. [25| and experimentally studied in Refs. 



D. Figure-8 Curve 

In order to investigate the consequences of the non- 
Abelian structure of the Uhlmann holonomies of the L 
and S subsystems, we consider here the class of "figure-8" 
loops on the parameter sphere of magnetic field directions 
shown in Fig. [U for which the holonomies can be calcu- 
lated explicitly. These loops are chosen to enclose no net 
area. Since any Abelian geometric phase of a angular 
momentum system is proportional to the area enclosed 
on the parameter sphere any such phase must van- 
ish for the figure-8 loops; a fact that has been used to 
demonstrate the Abelian nature of the Berry phase ex- 
perimentally by using nuclear magnetic resonance tech- 
niques [27|. Similarly, the mixed-state geometric phases 
|13| and the Berry phases in the LS system are all zero 
due to their Abelian nature. In contrast, we show that 
the Uhlmann phases of the subsystems along this class of 
loops are in general non-zero, which is a clear signature 
of the non-Abelian structure of the underlying Uhlmann 
holonomy. 

One can parametrize the loop F = F" o F' in Fig. [T]as 



F'( 



{9'o{t) e [0,7r),</)o) 

{TT,(t>'o{t) e [(/)o,(/)i]) 

{e[{t) e (7r,O),0i) 
(O,0'i(Oe [01,01 +7r]) 



(27) 




FIG. 1: Figure-8 loop on the parameter sphere of magnetic 
field directions. The Abelian geometric phase of the whole 
system along this path is zero because the net enclosed area 
is zero. In contrast, the phases of the Wilson loop variable 
associated with the Uhlmann holonomy for the L and S are 
in general nonzero. 



and 



F"( 



T'l 

r" 

^ 2 

F" 
^ 3 



(0(,'(t)e [O,^),0i+7r) 

(^,0{,'(i)e [01+^,00 + 

(0,0'/ (t) G [0o+7r,0o]) 



(28) 



with {9Q{t), . . . , 01 (i)} being a time ordered set of smooth 
functions. 

For the extremal subspaces ^ = ifJ.c, the holonomies 
of the L, S, and J systems are trivial in the sense that 

'^uhi(C'x±)' ^ ~ ^t^tJ^ become projection operators. 
This follows from the Abelian nature of the extremal 
states and from the fact that the net area vanishes for 
F. 

On the other hand, in the case where |^| < I + \, 
the holonomies turn non-Abelian and the corresponding 
Uhlmann phases might be non-zero. We demonstrate 
this in detail for the L subsystem. 

By using that [/^(O, 0) = 1, integration along F' and 
F" yields 



-bl al 

al bL 
~bl aL 



(29) 



Here, 



aL = e" 
bL 



=^W-^o)cos2(^)+sin2(^) 
where xl = wC^^'^^i:. Consequently 



\bL 



(30) 



(31) 
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where C 



[1,1^) 

L,± 



'■L.,± 



We may interpret the Uhlmann holonomy in Eq. pip 
as a rotation in an abstract space defined by the two 
states |/,/i± 5). Exphcitly, if we let = tanle"'", 
and introduce the effective Pauli operators ct£ = + 
— ^1 + h.c, a\ = —i\l,fJ, + ^){l,fJ- ^ 5I + h.c, and 
f^L = ^1 - 5I, defining an 

internal xyz coordinate system, the Uhlmann holonomy 
in Eq. (I3ip can be viewed as a rotation by an angle k 
around an axis in the xy plane making an angle rj with 
the X axis in the internal space. 

By applying the definition in Eq. ([Sj to the Uhlmann 
holonomy given by Eq. pip , we obtain the Uhlmann 
phase as 



IL;± 



arg 



Tr [/, 



^uhi(<f)) 



argC, 



where 



^ = 2[cos ((^1 - (^0) sin^ (xl) + cos^ (xl)]- 



(32) 



(33) 



Hence, 



loops for which the Uhlmann phase j^!^ is tt, in contrast 
to the corresponding mixed-state geometric phases that 
always vanish for such loops. This result is due to the 
non-Abelian nature of the underlying Uhlmann holonomy 
for non-extremal states. The existence of non-zero J^l!^^ 
is furthermore related to entanglement in a non-trivial 
way: one can show that a non-zero 7^^'.'^'' may occur only 



if 



< 



< 



One may verify that the phase 75.'^'' of the Wilson loop 
variable associated with the Uhlmann holonomy of the S 
subsystem may similarly be tt for certain figure-8 loops, 
while or undefined for other loops. Note in particular 
that the necessary condition on concurrence for a non- 
zero Jg!^^ now reads i < C^'''^^ < |, which is different 
from the above condition for the L subsystem. Thus, 
there may be energy eigenstates that allows for a non- 
zero 75.'^^^ while 7^'.'^'' must be zero and viceversa. In fact, 

if (Z + ■^)^ — /i^ > 9, then only one of the two Uhlmann 
phases can be tt for any given loop on the parameter 
sphere. 



7i'f = <; c < 



0, ^>0 

undefined, ^ = 0. 



(34) 



The points in the space (0i — 0o,xl oc C^-^^) where 
^ vanishes form a nodal line along which the Uhlmann 
phase 7^^'.'^^ is undefined. The points along this line are 
analogous to the nodal points found in Ref. |ll| in the 
case of the mixed-state geometric phase for this system. 




*,-*o 

FIG. 2: (Color online) ^ defined in Eq. (I33p as a function 
of — 4>o, for I = S, jj. — |, and tfie coupling strengths 

(? = 3, 13, 37, 50. The corresponding phase 7£.'^'' of the Wilson 
loop variable associated with the Uhlmann holonomy of the 
L subsystem is given arg^. 

In Fig. [2] we have plotted ^, as a function of (f>i — (f>o for 
/ = 3, ^ = |, and g = 3, 13, 37, 50. The figure shows that 
^ can in fact satisfy each of the three possible conditions 
displayed in Eq. (IMl) . In other words, there exist figure-8 



E. Additivity 

In this section, we explore the relation between the 
Uhlmann phases of the whole system and the subsystems. 
We restrict to cyclic evolutions, for which the Uhlmann 
holonomy takes the Wilson loop form. 

For fi — ±Hc, Eq- (1^ yields the sum rule 7j.^° = 

7^.^'^ +7s. J? for the corresponding Uhlmann phases. This 
confirms the expected sum rule for the geometric phases 
of product states. 



For \fi\ < I + ^, we put A^l'^'^ = A 



where A 



7j;± 



Cm) 

L;U(i; 



i;U(i) 



A 



L;SU(2) 



= /i(l — cos6)d4>, and obtain 



7L;± 



-arg 

7L;± 



" 7S;± 

Tr {Ve 



Z-;SU(2) 



)-( 



Pe 



A7 



(35) 



Since the trace of an SU(2) matrix is real, the deviation 



■ 1^1'. ± + 7s- ±^ can only 



from the sum rule 7j.'^'* 
be TT for a cyclic evolution in this model system. Similar 
to the nodal points found in Ref. for the mixed 

state geometric phases [l^ of the L and S subsystems 

there exist loops for which either Tr ^Pe 
Tr (^Pe"'^r'4L';su(2)^ vanish so that A7('''') becomes un- 
defined. These loops are nodal points of 7^^'.'^^'' or 75.'^'', 
respectively. 

We demonstrate that A7^'^^'' can be zero, non-zero, or 
undefined for the "orange slice" loop defined in Eq. (|27p , 
which connects the two poles on the parameter sphere 
twice along geodesies at (j)^ and 4>i . 
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The holonomy of the L subsystem for the orange shce 
loop r' can be found in Eq. (P^ . For the S subsystem, 
we obtain 



a*g -b*s 
bs as 



(36) 



where as and 65 are obtained from and 6^ in Eq. ([50)1 
by replacing xl with xs = C*^''^V. The corresponding 
Uhlmann phases read 

7^^ = -A*2(<^i - 0o), 
j^l:^^ = -^^2(01-00) 
+ arg 



cos2 (2^)cos((/.i-(/.o) + sin2(^) 



arg 



cos2 (^) cos ((/.I - 0o) + sin2 (^) 



(37) 



Thus, A^^'''^^ = — argC, where 



c 



cos2(^)cos(0i-0o)+sin2 (^) 



cos^ (^) cos (01 



))+sin^ (^)].(38) 



is if C > 0, 77 if C < 0, and undefined if C = 0. 
All these three cases are visible in Fig. ([3]), in which we 
have plotted as a function of (pi — 0o for fixed I, fi, 
and 5 = 3, 20, 50. Figure [3] confirms that there are some 
loops on the parameter sphere for which A7''^''' 7^ 0, and 
therefore tJ'^^ ^ j^^!^^ + 75;± • 

To compare with the mixed state geometric phase pro- 
posed in Ref. [l^, we first consider the spectral de- 
composition {pk,\tpk)} of an arbitrary non-degenerate 
density operator p that undergoes cyclic unitary evolu- 
tion. The resulting mixed state geometric phase /3 reads 
/? = arg (^J2k Pke^^'°), where Pk is the cycfic pure state ge- 
ometric phase of {ipk)- Using this expression, the mixed- 
state geometric phases I3(C^'^) and (3{Cg'!^) of the L 
and S subsystems read [ill 



^(C[' ±^) = ± arctan ( cosa^'^^^ tan - ) , 



=F arctan I cosa*-''^-* tan 



which implies the sum rule /3(C[';^^)-H/3(C^''^^) = -^n = 

P{Cj'^^) = 7j.±''. Thus, the sum rule for /3 is satisfied for 
any loop F on the parameter sphere. This again demon- 
strates the difference between the phase of the Wilson 
loop variable associated with the Uhlmann holon omy and 
the mixed-state geometric phase defined in Ref. [L^ . 

We note that, the alternative Uhlmann phase tpuhi de- 
fined in Eq. |4]is in general non-zero for the figure-8 path 
in Fig. [TJ as a consequence of the non-Abelian nature of 
the underlying Uhlmann holonomy. On the other hand, 
the corresponding deviation Aip^^^^^ from the sum rule 



(39) 



0.5 



-0.5 




-g=3 

-g=20 

g=50 



1 



FIG. 3: (Color online) defined in Eq. (|38} as a function 
of 01 - 00 for / = 2, = i, and g = 3,20,50. The dif- 
ference A7(''f' = 7^';^^ - 7['.'^> - 7^';^' of the phases of the 
corresponding Wilson loop variables is given by — arg (. 



is no longer restricted to or tt, but can take any real 
value for cyclic evolutions. Thus, in general there is a 
difference between i^Suhi and 7, and they both differ, in 
the sense of the sum rule, from the mixed-state geometric 
phase in Ref. [l3j . 



V. CONCLUSIONS 

Uhlmann's quantum holonomy along density operators 
is a concept that allows for studies of geometric phases 
of general quantum states undergoing arbitrary quantum 
evolutions. Its relevance to various aspects of physics 
have been demonstrated in the pastj_such as for open sys- 
tem evolutio n 12811 , interferometry [29] , many- body quan- 
tum systems [30], as well as for Yang- Mills theory 31) and 
Thomas precession in special relativity [12] • These works 
have been triggered in part by the need to test the conjec- 
tured resilience of holonomic quantum gates for quantum 
computation to various kinds of errors, such as noise and 
decoherence. Recently, a first experimental test of the 
related Uhlmann geometric phase, utilizing nuclear mag- 
netic resonance techniques, has been carried out (20j . 

Here, we have analyzed Uhlmann's quantum holon- 
omy by considering a physical model system in which 
the Uhlmann holonomies for LS'-coupled hydrogen-like 
atoms in a slowly rotating magnetic field have been 
computed. We have shown that the holonomy of the 
total angular momentum has Abelian structure. Fur- 
thermore, its corresponding phase is exactly the asso- 
ciated standard geometric phase [l6| for open or closed 
paths on the parameter sphere of magnetic field direc- 
tions. For the holonomies of the L and S subsystems, 
we have shown that, in analogy with the Levay geomet- 
ric phase defined for two-qubit systems [2^, depending 
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on whether the energy eigenstate of the whole system is 
a product state or an entangled state, the correspond- 
ing holonomies are Abelian or non-Abelian, respectively. 
In the case of entangled states, there is an explicit de- 
pendence o^the gauge field vector potential on the con- 
currence [2^, which interpolates the standard Abelian 
(Berry) and non-Abelian (Wilczek-Zee) cases. In other 
words, our analysis demonstrates that the rich geomet- 
rical nature of the Uhlmann holonomy incorporates as a 
limiting case the Wilczek-Zee holonomy, which is char- 
acterized by maximum quantum entanglement between 
the L and S subsystems. 

In the analysis of the phase of the Wilson loop vari- 
able associated with the Uhlmann holonomy, we have 
pointed out that this phase, unlike the mixed-state geo- 
metric phase p^ . possesses a non-Abelian structure and 
may therefore be non-zero even for loops on the param- 
eter sphere that enclose no net area. We have also elu- 
cidated that the phases of the Wilson loop variables of 
the corresponding Uhlmann holonomies for the L and S 
subsystems add up to that of the whole system for spe- 
cific paths; for other paths the sum may differ by tt from 
the Berry phase of the whole system. 

Furthermore, we would like to point out that previ- 
ous theoretical [T7l - [l9| and experimental [13] work ana- 
lyzing the relation between the Uhlmann holonomy and 
the mixed-state geometric phase all employ a notion of 
Uhlmann phase [see Eq. (|4])] which differs both concep- 
tually and numerically from the phase of the Wilson loop 
variable associated with the Uhlmann holonomy, as de- 
fined in Eq. ([S]). This alternative concept is the phase 
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